Quantum Transport of a ID Degenerate Bose Gas in a Lattice: The Role of Different 

Quantum Fluctuations 
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We study the effect of quantum fluctuations on the quantum dynamics of ID Bose gases in 
an optical lattice using different basis expansions and the Truncated Wigner approximation. The 
initial spatial distributions of the quantum fluctuations are very different when using a limited 
number of plane- wave (PW), single-harmonic-oscillator (SHO), and the self-consistently determined 
Bogoliubov (SCB) modes. The transport properties of the Bose gas, characterized by the phase 
coherence and number fluctuations, in the PW basis are distinct from those gained using SHO 
and the SCB basis. The calculations using the SCB modes predict greater phase decoherence and 
stronger number fluctuations than the other choices. Furthermore, we observe that the use of PW 
modes overestimates the extent to which atoms are expelled from the core of the cloud, while the 
use of the other modes only breaks the cloud structure slightly which is in agreement with the 
experimental observations [1[. 

PACS numbers: 03.75.Dg, 03.75.Lm, 39.20.+q 



I. INTRODUCTION 

The influence of quantum fluctuations on ultracold 
atomic gases have attracted much attention in recent 
years. In periodic optical lattice (OL), bosons can un- 
dergo a superfluid-to-Mott-insulator (SF-MI) quantum 
phase transition because of the loss of phase coherence 
induced by quantum fluctuations if the tunneling became 
small enough 0, E| . Theoretical work has shown that for 
a system where condensate oscillations are overdampcd 
in a ID OL, there is a smooth crossover between the 
classical transition and the SF-MI phase transition at 
zero displacement as the quantum fluctuations increase 
In addition, the damped motion of a ID Bose gas in 
a shallow OL has been demonstrated in experiment [l[. 
Many theoretical studies have indicated that quantum 
fluctuations can result in strong dissipation which leads 
to dynamical instabilities [3-0] . 

The truncated Wigner approximation (TWA) devel- 
oped in the field of quantum optics provides an attrac- 
tive framework to explore the effect of the quantum fluc- 
tuations on the dynamics of Bose-Einstein condensates 
(BECs) [|, [Ml. The fundamental idea of this method 
is to expand quantum field operators in the equation of 
motion for a Bose gas system in the Wigner represen- 
tation and, by the neglect of the third-order functional 
derivatives in the equation, derive a generalized Fokker- 
Planck equation (FPE) for a classical field problem. This 
results in a stochastic differential equation for the clas- 
sical field ^(x, f), which coincides with the mean field 
Gross-Pitaevskii equation [j| [TJ| 

ififittffot) = i/ *(x,i)+ 5 |*(x,i)| 2 *(x,t). (1) 

where H = -h 2 d^/2m + V(x,t). The quantum- 
mechanical vacuum fluctuations are simulated by adding 
appropriate classical fluctuations in addition to the co- 



herent field of the initial state of the BEC. Specifi- 
cally, the quantum fluctuations are introduced in the ini- 
tial condition for the classical field \&(x,i), by the sum 
^(x, t = 0) = V( x ) + £( x ); where ip(x) and £(x) are 
the real and virtual particle fields, respectively. The real 
particle field is the suitably chosen initial form for the 
classical field. The field of virtual particles is decom- 
posed into a basis set £(x) = Ylj=i £j0j( x )j where 4>{x) 
is the mode function and the amplitude in each mode is 
Gaussian with the properties, = ^i.j, (CiCi) = 0- 

Expectation values are evaluated by averaging over an 
ensemble of mean-field calculations with random initial 
conditions. It implies that the stochastic element only 
appears in the initial conditions. 

Recently, the TWA has predicted quantum turbulence 
in the collisions of two condensates [13[ and explained 
the disruption of reflecting BECs from Silicon surface 
[l3 |. In these theoretical works, zero-temperature quan- 
tum fluctuations are modeled by adding, on average, a 
half quantum into each low-energy PW mode. There- 
fore, the quantum fluctuations added into the nonuni- 
form condensates are spatially uniform, suggesting that 
the distribution of quantum fluctuations docs not depend 
on the geometry of atom clouds and interatomic interac- 
tion. However, if the quantum is added into a limited 
number of distinct modes in the treatment of TWA, the 
corresponding initial wave functions of the classical field 
are different. In addition, the wave vector is not a good 
quantum number for a harmonically trapped condensates 
Oil 03 1 except for a ID, very elongated, condensate fl5j . 
In the frame of Bogoliubov theory, the geometry and in- 
teratomic interaction of condensate clouds have a crit- 
ical effect on the distribution of quantum fluctuations. 
The TWA with Bogoliubov theory in the treatment of 
quantum fluctuations has shown reasonable results in the 
nonadiabatic dynamics and dissipative dynamics of quan- 
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turn transport in optical lattices Q . Provided a complete 
set of modes are used, the choice of basis would be im- 
material from the point of view of the validity of the 
method, and the computationally simplest basis is pre- 
ferred. However, the method requires a cut-off to be im- 
posed on the mode expansion which distinguishes "low 
energy" modes, that properly contribute to fluctuations 
in the condensate, from "high energy" modes that do not. 
The use of a truncated basis means that it is no longer 
clear that all basis sets are equivalent. 

In this paper, we consider three choices of basis set: 
PW modes, SHO eigenstates, and SCB modes, to con- 
struct the quantum fluctuations with different initial 
states in TWA. Our purpose is to examine the effect of 
these different quantum fluctuations on the dynamics of 
a strongly confined Bose gas, oscillating in a very shal- 
low ID OL. While we find, as shown in Ref. (HQ, that 
the inclusion of quantum fluctuations does indeed lead to 
damping of the dipole oscillations of the center-of-mass 
(cm.) of the condensate, we note that the choice of re- 
stricted basis set has a qualitative effect on the evolution 
of the spatial density distribution of the cloud, the loss 
of phase coherence, and the number fluctuations. These 
qualitative differences originate from the initial distribu- 
tion of the quantum fluctuations. In the PW basis, some 
atoms are kicked out from the core part of atomic cloud, 
resulting in a density distribution with a long tail. This 
does not agree completely with the experimental observa- 
tions [l| , where the cloud in damped transport has a sim- 
ilar width to the undamped case. However, there is only 
slight disruption of the cloud structure in the other two 
cases. Moreover, by using the same number of modes, 
Bogoliubov theory offers stronger damping of the cm. 
trajectory than others. 

The remainder of this paper is organized as follows. In 
section II we summarize the truncated Wigner scheme 
and our specific implementation of it for the three choices 
of basis states. In section III we describe the results of 
the simulations and highlight the differences in dynamical 
evolution between the different choices of basis set. In 
section IV we summarize these and draw more general 
conclusions about such simulations. 



II. THEORETIC MODEL AND NUMERICAL 
METHODS 

Our theoretical model is based on the recent experi- 
ment [l|, where damped dipole oscillations of a ID Bose 
gas in a shallow OL were observed. In this experiment, 
the ID Bose gas is formed in an array of independent 
"tubes" , produced by applying a strong transverse 2D 
OL potential to confine a trapped 3D atomic condensate. 
Then the tubes are corrugated adiabatically by using a 
very shallow ID lattice along the axial direction. The 
dipole oscillations of atoms along the weak axial lattice 
were excited by suddenly displacing the harmonic trap 
and the cm. velocity was imaged. 



The transport dynamics of the system based on the 
above model is studied by using the TWA method. In a 
very shallow OL, the dynamics of the system can be de- 
scribed by the ID GPE (Eq.([T])) with the coupling con- 
stant g = g\D — 2huj±a. We use the s-wave scattering 
length a = 5nm for 87 i?6 and us± — 2ir x 38kHz, origi- 
nated from the quadratic expansion of the strong trans- 
verse 2D OL around the local minima [H,i]. The poten- 
tial energy profile of the ID optical lattice in the axial 
direction is characterized by VoL(x,t) = A(t) sin 2 (wx / d) , 
where d = 405nm with respect to the laser wavelength 
A = 810nm. The amplitude A(t) is assumed to zero 
initially, ramped gradually up to E r /2 in 2.65ms as 
exp(fci) — 1 where E r — h 2 /2m\ 2 is the photon recoil 
energy, and k is an constant determined from the ramp- 
ing time. Thereafter A(t) remains unchanged and the to- 
tal confining potential is V(x, t) — Vol(x, t) + muj 2 x 2 /2, 
where u) = 2tt x 60Hz is angular frequency of the har- 
monic trap. At time t = 2.65ms, we abruptly displace 
the harmonic trap through a distance Ax — 3/im, accel- 
erating the Bose gas in the OL. 

We consider the dynamics of the system for distinct ini- 
tial states by introducing different mode functions <p{x) 
in the field of virtual particles £,(x). In case I, the con- 
densate is assumed to propagate in the OL in absence of 
quantum fluctuations, i.e. £*(x)t;(x) — 0. In case II, the 
initial amplitude of virtual particle field, is 1/2 

and (f>i( x ) = exp(ikix) / '\fL ', where ki is the wave vector 
in the i-th mode and L is the length of system [l3|) [lH ■ 
In case III, the mode function, ip n (x), is the eigenfunc- 
tion of the SHO, resulting in the quantum fluctuations 
in a limited number of energy modes. In case IV, we use 
Bogoliubov approximation 

$(x,t) = ip (x)a + ^2[u j (x)a j (t) - v](x)d](t)}, ( 2 ) 

where ipo(x) is the ground state solution of the GPE, the 
number of atoms < ajao >= No, and aj, dj are quasipar- 
ticle field operators. Then we diagonalize the Bogoliubov 
equations, 

Cuj (x) - gi D N ^ 2 {xjvj (x) = E 3 u (x) 
Cvj{x) - giDN ^Q 2 (x)uj(x) = -EjVj(x), 

in the subspace orthogonal to the condensate space to 
obtain the quasiparticle amplitudes, Uj{x), Vj(x) and ex- 
citation energy Ej. In Eqs.©, C = —h 2 d 2 /2m + U e ff{x) 
with 

U Bff (x) = V(x) + 2g 1D N Q \Mx)\ 2 ~ V (4) 

the effective potential of quasiparticles [l8j]. The initial 
state of fy(x) can be obtained using Eq.Q and dj, d^ 
(J > 0) are replaced by the random variables (ay, a*), 
obtained by sampling the corresponding Wigner distri- 
bution of the quasiparticles in a thermal bath [l9j , 

2 

W(ctj,a*) — —tanh(i]j)exp[—2\aj\ 2 tanh(rij)], (5) 
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FIG. 1: The spatial distribution of density fluctuations, at- 
tributed respectively by adding a half quantum into each PW 
(a), SHO (b), and SCB mode (c). The effective potential U e ff 
defined by Eq.gJ) (d). 
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where rjj = Ej/2ksT. At temperature T — 0, the quan- 
tum fluctuations are described by t;{x) = J2j u j( x ) a j ~ 
v 3 (x)a* in the SCB modes based on Eq. © andEq.©. In 
this paper, we consider the total atom number iVo = 100, 
the number of modes, M < 80, and always restrict our 
discussions at zero temperature where the Bogoliubov 
approximation can be justified. Two groundstates with 
optical lattice amplitudes A = and E r /2 are calculated 
numerically by evolving the GPE in imaginary time [2fj| . 

The plane waves are the simplest computationally as 
the field of virtual particles £(x) is obtained using the 
Fast Fourier transform method. The SHO basis has the 
advantage that the support of each mode is weighted 
towards the region in which the cold atoms are likely 
to be found and, while the construction of £(x) is more 
complicated, the modes are at least pre-defined and are 
relatively straightforward to generate iteratively. The 
SCB modes should be good approximations to the quasi- 
particle excitations of the system, at least for weak inter- 
actions, and are well adapted to the specifics of the po- 
tential in which the cloud moves. The calculation of £(x) 
in more computationally demanding because the form of 
the modes has to be determined by diagonalization of a 
large matrix whose elements depend on the form chosen 
for the initial order parameter field, ipo(x). 

To determine the validity and precision of the results, 
different numerical methods are applied in our simula- 
tion. In case I, we obtain dynamic results by using 
the Crank-Nicolson method to solve the time-dependent 
Gross-Pitaevskii equation. In case II, a distinctive nu- 
merical method, RK4IP-P [H, [HI, is utilized. Both 
methods are used in case III-IV to examine their equiv- 
alence. Solutions of the cm. trajectory, phase coherence 
and number fluctuation in TWA, are calculated numer- 
ically using different numbers of realizations. We find 
that some results can be significantly different when the 
number of realizations is small, i.e. less than 50. There- 
fore, we show all results based on 200 realizations for 



FIG. 2: (Color online) The time evolution of the cm. motion 
of the Bose gas in cases I (solid black line), II (solid green line), 
III (dashed red line), and IV (solid blue line). The fitting of 
the cm. trajectories in cases III and IV (dotted black line 
in (b) and (d)) from Eq.©. The corresponding parameters: 
B = 0.22^m, 0.1/xm and 7 = 22.89,10.44 for cases III and 
IV, respectively , while Q — 2tv x 60 and A = 3.0/im for both 
cases. 



numerical consistency. 



III. RESULTS AND DISCUSSION 

In FiglU we show the initial spatial distributions of the 
quantum fluctuations (£* (x)^(x)} with respect to three 
sets of basis modes. In the PW modes, (£*(x)£(x)) = 

JZ Sfj=i (£i*£j) ex P[*(^i — kj) x ]> indicating a uniform spa- 
tial distribution of the quantum fluctuations, while the 
limited number of low-energy SHO modes are dominated 
by the Gaussian function, which induces a localized spa- 
tial distribution of the quantum fluctuations. In contrast, 
the quantum fluctuations in the SCB modes give a spatial 
distribution with a double peak beside the atom cloud, 
which coincides with the effective potential of the quasi- 
particles Eq.Q [lfl, as shown in FigHJd). 

Qualitative insight into the properties of the quantum 
dynamics can be gained by using classical dissipative dy- 
namics to compare with our quantum simulations. We 
model the cm. motion as a damped harmonic oscillator 

m 
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where the cm. displacement is defined by 

f^x^ix^dx 



X c . m . = (x)(t) = 



fZoM^dx ' 
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FIG. 3: (Color online) Phase coherence C3 in different cases 
(a), (b), and Cj in case III (curves from top C3, C5, Cg, C10) 
(c). 
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FIG. 4: (Color online) Number fluctuation Ano in different 
cases (a), (b), and An,- in case III (curves from top Ano, 
Anio, Ari2o, Aruo) (c). 



In underdamped case, Eq.fJB]) have a solution 

,- 7 (t-to) [AcosQ,(t - t ) + Bsmfl(t - t )} 



(8) 



7 , and to determined by 



with B = jA/Q, tt 
the initial phase. 

In Fig|2j we show simulation results of the cm. trajec- 
tories for case I- IV in comparison with Eq.©. In case I, 
there is undamped motion of the cm. while the visibly 
damped trajectories occur in presence of the quantum 
fluctuations in the other cases. The damped trajectories 
(x)(t) cannot be fitted by Eq.©, especially at long times, 
highlighting the complexity of the effect of the quantum 
fluctuations on the transport of the system. Therefore, 
the chosen fitting parameters satisfy the numerical calcu- 
lations well in the first period but hardly characterize the 
subsequent behavior of the damped motion. To compare 
the time-dependent damping rate in different cases, we 
introduce A = \x\{D$/ D\) [4|, where Dq and Di are the 
cm. positions at t — 10.9ms and 27.6ms ( FigJSJc)). Our 
calculations for cases III and IV (Fig|2^c) ) give A = 0.21 
and 0.26 respectively, therefore indicating that, with the 
same number of basis states, the quantum fluctuations in 
SCB modes inhibit more strongly the motion of the Bose 
gas than in other modes. There are slight differences be- 
tween the cm. trajectories in cases II and III, as shown 
in FiglUJa). 



The damped motion of the cm. suggests that there 
might be a loss of phase coherence, and number fluctua- 
tions in the system. In order to avoid the complications 
arising from the symmetrically ordered multimodc field 
VP |5|, we define the ground state operators bj for each 
individual lattice site j: 



b 3 (t) 



j th well 



(9) 



where ^(x, t) is the stochastic field, determined by Eq.([T]) 
and ^o{x) is the groundstate wave function in the com- 
bined harmonic trap and OL. The normalized phase 
coherence between the central well and its j'th neigh- 
bor and the atom number fluctuations in the jth site 
are separately set as Cj — | ( Sq £>j- ) | / y/norij and Arij = 

In Fig|3]and FigUl we show respectively the configura- 
tions of the dynamic phase coherence and number fluc- 
tuations for cases II-IV ((a) and (b)) and on different 
spatial sites for case III ((c)). In case III, there is an 
increase of the phase coherence C3 at t < 4.8ms and 
then its loss at t > 4.8ms in the first period of the cm. 
oscillation, while the number fluctuations A713 decrease 
and grow correspondingly, displaying the opposite dy- 
namic behavior between the phase coherence and num- 
ber fluctuations at nearby sites. This behavior demon- 
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FIG. 5: Density distribution at t = 8.65ms in cases II (a), III 
(b) and IV (c). It — yjh/mu is the length unit. 



strates that the quantum fluctuations tend to suppress 
the phase coherence. Moreover, the phase coherence Cj 
and number fluctuation Arij are correlated essentially to 
the damped trajectory of the cm.. At t = 2.65ms, the 
cm. in Fig [2] is far from the center of the trap potential, 
and thus C3 and An , as shown in Fig(3f c) and FigSJc), 
identify the properties in the outer region of the Bose 
gas where the number fluctuation is large and the phase 
coherence is weak. Then the Bose gas moves toward the 
center of the trap, and Auq decreases gradually while C3 
increases. 

However, qualitative differences are particularly 
marked at short times, as depicted in FigJ3Ja) and 
FigHJa), indicating at least a marked dependence of the 
"equilibration time" for such simulations on the choice 
of basis set. This is most marked for the PW basis in 
which the short-time oscillations in the fluctuation statis- 
tics are quite different to the others and only approach 
them at times rather longer than the period of the dipole 
oscillations. FigO (b) and Fig0] (b) show that the short- 
time variation of C3 and Ano based on the SCB modes 
is more complicated than the other two modes. This 
might originate from the greater complexity of the ini- 
tial distribution of the quantum fluctuations in the SCB 
modes. Comparing the phase coherence for cases III and 
IV, shown in Figf^b), we find that the time average of 
C 3 from 2.65ms to 42.65ms for case IV is 0.63, 5.21% 
smaller than that for case III (0.66). The corresponding 
mean value of the number fluctuation in the time inter- 
val for case IV is 3.95, a little higher than 3.88 for case 
III. These indicate that with the same number of basis 



modes, the quantum fluctuations in the SCB modes cause 
stronger phase decoherence and number fluctuation than 
those in the PW and SHO modes. 

Interestingly, we also find that, during the transport 
process, some atoms are kicked out from the core re- 
gion of the atom cloud in case II (FigfSJa)) while the 
quantum fluctuations only break the inner configuration 
of the condensate slightly in cases III (Fig(5jb)) and IV 
(FiglSJc)). In the experiment [l|, they did not observe a 
significant difference in the time of flight width between 
atoms that undergo damped harmonic motion and those 
that are unexcited but held for an equal time. 

Of course, we can strictly only say that the three basis 
sets produce different results, but it is natural to assume 
that the simulations using the SCB modes have more in- 
formation about the specific physics of the system built 
into them and should be the most reliable. We then 
note that for some properties, the gain in using the SCB 
modes compared to the SHO modes is small given the 
considerable increase in computational complexity. How- 
ever we note that in the simulation that uses PW modes 
for trapped BECs, a long equilibration time should be 
allowed before any non-trivial dynamical processes are 
allowed to occur in the simulation to avoid spurious tran- 
sient effects. 



IV. CONCLUSION 

In conclusion, we have studied the effect of different 
choices of basis set for the inclusion of quantum fluc- 
tuations on a ID Bose gas in an OL using the TWA. 
Specifically we have used PW, SHO and the SCB modes 
to decompose the quantum fluctuations and examine the 
dynamics of the system. The difference in the predic- 
tions of the phase coherence, number fluctuations and 
density variations at short time indicates that the choice 
of the basis set has a substantial influence on the qualita- 
tive features of the transport. The use of the SCB modes, 
which gives greater phase decoherence and stronger num- 
ber fluctuations on the damping dynamics than the other 
two choices, incorporates more of the underlying physics 
of the problem and presumably leads to the most reli- 
able results. For the trapped condensate, the SHO basis 
give results that are qualitatively similar to those from 
the SCB modes, at rather less computational effort. It's 
therefore clear that it is important to select the basis set 
that best suits the particular system to be studied. 
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